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Abstract - Analytic expressions for the magnetic induction and 
its spatial derivatives for a circular loop carrying a static current 
are presented in Cartesian, spherical and cylindrical coordinates. 
The solutions are exact throughout all space outside the 
conductor. 

Index Terms - Circular Current Loop, Magnetic Field, Spatial 
Derivatives. 

I. Introduction 


as shown (Fig. L). It is assumed that the cross section of the 
conductor is negligible. 



A nalytic expressions for the magnetic induction (magnetic 
flux density, B) of a simple planar circular current loop 
have been published in Cartesian and cylindrical coordinates 
[1,2], and are also known implicitly in spherical coordinates 
[3], In this paper, we present explicit analytic expressions for 
B and its spatial derivatives in Cartesian, cylindrical and 
spherical coordinates for a filamentary current loop. These 
results were obtained with extensive use of Mathematical 
and are exact throughout all space outside of the conductor. 
The field expressions reduce to the well-known limiting cases 
and satisfy V • B = 0 and V x B = 0 outside the conductor. 

These results are general and applicable to any model using 
filamentary circular current loops. Solenoids of arbitrary size 
may be easily modeled by approximating the total magnetic 
induction as the sum of those for the individual loops [4]. The 
inclusion of the spatial derivatives expands their utility to 
magnetohydrodynamics where the derivatives are required. 

The equations can be coded into any high-level 
programming language. It is necessary to numerically evaluate 
complete elliptic integrals of the first and second kind, but this 
capability is now available with most programming packages. 


Fig. 1. Circular current loop geometry. 

The vector potential , A, of the loop is given by [3]: 
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where r , 6 , and (p are the usual spherical coordinates, and the 
argument of the elliptic integrals is 
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Note that we use k 2 for the argument of the elliptic integrals. 
This choice is consistent with the convention of Abramowitz 
and Stegun [5] where m = k 2 . 

For a static field with constant current, the B components in 
spherical coordinates are [3]: 


II. Spherical Coordinates 

We start with spherical coordinates because this is the 
system usually used in the standard texts. The Cartesian and 
cylindrical results in Sections III and IV were derived from the 
spherical coordinate results. 

The current loop has radius a , is located in the jc-y plane, 
centered at the origin, and carries a current / which is positive 
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Analytic expression for the field components and their 
derivatives in spherical coordinates are given below. For 
simplicity we use the following substitutions: 
a 2 =a 2 + r 2 -2ars\n & , f$ 2 = a 1 + r 2 3- 2arsin 6 ,k 2 = l -a 2 /p 2 , and 
C = p {) I/tt . 
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Field Components: 


Field Components: 
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Spatial Derivatives of the Field Components: 
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III. Cartesian Coordinates 

The field components and their derivatives in Cartesian 
coordinates are given below. These are easier to use when 
rotations or translations are needed and obviate the need to 
transform the basis vectors. The following substitutions are 
used for simplicity: p 2 = x 2 + y 2 , r 2 = x 2 + y 2 +z 2 , 

a - = cr + -2cip , p~ = a" + r" + lap < k = 1 — a /p , y = jc -y , 
and C s p {) 1 /k . Note that p and r are non-negative. 
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Spatial Derivatives of the Field 
Components: 
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We have presented simple, closed-form algebraic formulas 
for the magnetic induction and its spatial derivatives of a 
filamentary current loop that are exact everywhere in space 
outside the conductor. Although these formulas are exact, they 
do require the numerical evaluation of elliptic integrals. 

Solenoids with circular cross sections of arbitrary size and 
configuration can be modeled by simply summing the 
contributions of each individual loop. 

There are, of course, other ways to obtain B for the basic 
circular current loop. For example, series expansions are 
available [3] and numerical integration via a finite element 
approach can be performed [6]. However, these suffer from 
limitations such as truncating the series expansions after some 
tolerance is reached or accepting some graininess when using a 
discrete grid. Our approach has neither of these limitations and 
yields results are that exact up to the limitations of the 
numerical arithmetic and the elliptic integral routines. 

The inclusion of the spatial derivatives allows convective 
derivatives to be found exactly and may prove useful for 
magnetohydrodynamics applications. 
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